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Abstract 

We study a class of piecewise linear solutions to the inviscid Burgers equa- 
tion driven by a linear forcing term. Inspired by the analogy with peakons, 
we think of these solutions as being made up of solitons situated at the 
breakpoints. We derive and solve ODEs governing the soliton dynam- 
ics, first for continuous solutions, and then for more general shock wave 
solutions with discontinuities. We show that triple collisions of solitons 
cannot take place for continuous solutions, but give an example of a triple 
collision in the presence of a shock. 

1 Introduction 

The subject of this paper is piecewise linear solutions of the PDE 

(u t + uu x ) xx = 0, (1.1) 

which we earlier [T] have called the derivative Burgers equation. This name 
refers of course to the well-known Burgers equation u t -I- uu x — vu xx and its 
special case the inviscid Burgers equation u t + uu x = 0, which is the prototype 
equation for studying shock wave solutions of hyperbolic conservation laws. In 
some applications one considers also forced Burgers equations with terms of 
the form F(x,t) on the right-hand side, often written as F = —dV/dx with a 
potential V . Since equation (jl.ip is equivalent to u t + uu x = A{t)x + B(t), it is 
perhaps more appropriate to talk about it as a forced inviscid Burgers equation 
with linear force (or quadratic potential). Moreover, the latter equation can be 
rewritten as 

ut + \{u 2 ) x = A(t)x + B(t), (1.2) 
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which makes sense for a much larger class of functions than just u £ C 1 (R). 
For example, if u € Lf oc (R) we can interpret (|1.2p to hold in the sense of dis- 
tributions. One could work with distributions acting on test functions tp(x, t) £ 
£>(R 2 ), but the following simpler interpretation is sufficient for our purposes 
here: we view u(x, i) as a mapping that takes a real number t to a function 
u(-,t) £ L 2 oc (R) which we can identify with a distribution in £>'(R). The 
derivative with respect to x is then the distributional derivative defined by its 
action on a test function if)(x) £ 2?(R) in the usual way, (u x ,ip) = — {u,ip x ), 
while the derivative with respect to t is the limit of a difference quotient. If 
equation p.2[) is satisfied in X>'(R) for each t, then we then say that it holds in 
a weak sense and that u is its weak solution. 

We were led to the Burgers equation by our previous work on peakon and 
shockpeakon solution of the Degasperis-Procesi (DP) equation 

u t - u txx + Auu x = 3u x ii xx + uu xxx , (1.3) 

an integrable wave equation discovered a few years ago [3J. Indeed, the 
problems treated in this paper are to some extent "toy problems" , but we hope 
that they might provide some guidance and intuition for the future study of the 
DP equation. 

Equation (|1.1|) can be obtained formally from the DP equation by substi- 
tuting x i— > ex, t t— » st, and then letting e — > 0. This "high-frequency limit" is 
a natural thing to try on the DP equation, since it is the same procedure that 
takes the celebrated integrable Camassa-Holm (CH) shallow water equation 0], 

ut — u txx + 3uu x = 2u x u xx + uu xxx , (1.4) 

to the Hunter-Saxton (HS) equation for nematic liquid crystals [5j[6], 

(u t + uu x ) xx = u x u xx . (1.5) 

The CH and DP equations both admit peakon solutions, which are multi- 
soliton solutions of the form 

n 

u(x,t) = 2jmfc(i) exp(- \x - x k (t)\), (1.6) 
fe=i 

where the functions 2Jfc(t) and m^it) (positions and momenta of the individual 
peak-shaped solitons) are required to satisfy a certain system of 2n ODEs in 
order for u(x, t) to satisfy the PDE in a weak sense. In shorthand notation 
these ODEs are Xk = u(xk), ihk = —(b — l)u x (xk), where b = 2 for the CH 
equation and b = 3 for the DP equation. One can think of this as an integrable 
mechanical system of n particles on the real line, simililar to, for example, the 
open Toda lattice. If follows from the rapid decay of e - ^' that iu = u(xk) ~ 
when all distances large, so it agrees with intuition to regard rrik as 

the momentum of the fcth particle. Asympotically (when t — > ±oo) the particles 
will spread apart, each moving with its own (nearly) constant velocity which is 
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nonzero and distinct from the other particles' velocities. The latter is a highly 
nontrivial fact for the DP equation pj Theorem 2.4]. 

There is an analogous class of solutions of the HS equation and the forced 
Burgers equation (|1.2p . namely the piecewise linear solutions 



In the shorthand notation used above, the governing ODEs take exactly the 
same form again: Xk = u(xk), fhk = — (6 — l)u x (xk), where 6 = 2 for the HS 
equation and 6 = 3 for the forced Burgers equation. However, for peakons the 
term rrike~ ][Xk ~ Xk \ usually dominates the other terms in the equation Xk = u(xk), 
while here we instead have the term rrik \x% — Xk\ which is zero while all other 
terms are large. Thus, in contrast to peakons where the interaction is strongly 
localized, these piecewise linear solitons influence each other more strongly the 
more separated they are. Although it is a bit hard to develop a useful intuition 
about these ODEs as a "mechanical" system (perhaps one can think of some 
kind of expanding gas with long-range correlations), the analogy with peakons 
still makes it natural to think of the piecewise linear solutions as being composed 
of some kind of solitons situated at the breakpoints Xk ■ (But we have not been 
able to make sense of the idea that the piecewise linear solutions are somehow 
high-frequency limits of peakons). 

In all four cases mentioned above, the ODEs governing the soliton dynamics 
can be explicitly solved using inverse spectral methods [3 [9l [3 [10l Q] [11] . In 
the forced Burgers case the ODEs are also easily solved directly by elementary 
methods, as we will see. 

In the Degasperis-Procesi equation (but not in the Camassa-Holm equation) 
there also appears a more complicated phenomenon, namely discontinuous so- 
lutions of the form 



u(x,t) = ^( TO fe (t) - s fe (i) sgn(x - x k {t)) J exp(- \x - x k (t)\). (1.8) 



Such shockpeakons [12] are governed by 3n ODEs for positions x k , momenta 
TOfc, and shock strengths s&. Even if one starts with the usual peakon ansatz 
(|1.6[) . shock solutions of the form (|1.8p can form after finite time when a peakon 
with mfe > collides with an antipeakon with mk+i < moving in the opposite 
direction. (In the CH equation, such collisions give rise to "zero-strength shocks" 
where u x momentarily blows up but u remains continuous, still being of the form 
(|1.6[) after the collision [8], and a similar thing occurs for the HS equation [6].) 
The shockpeakon ODEs have so far only been solved in the trivial case n = 1 
and in a very particular subcase when n = 2. The problem is that the Lax 
pair for the DP equation, which was crucial for deriving the peakon solution 
formulas, does not make sense for the weak formulation of the DP equation 
that is used when working with discontinuous solutions. 




(1.7) 



fc=i 
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The forced Burgers equation (|1.2p admits an analogous class of solutions, 
given by the discontinuous piecewise linear ansatz 

u(x,t) = TOfc(t) \x - x k (t)\ - s k (t) sgn(a; - x k (tj) J . (1.9) 
fe=i ^ ' 

Such solutions with shocks can form after finite time, even if the initial profile 
is continuous. Unlike the Degasperis-Procesi case, it turns out here that the 
extra generality of having jumps in u can be handled without problems. 

The outline of the paper is simple: we derive and solve the ODEs gover- 
ing piecewise linear solutions of the forced Burgers equation (| 1 . 2[> . first in the 
simpler case (|1.7|) of continuous solutions (using elementary methods and, for 
comparison, inverse spectral methods), then in the general case (11.91) of discon- 
tinuous solutions (by reduction to the previous case). We conclude with a few 
examples. 



2 Continuous piecewise linear solutions 

Theorem 1. The continuous piecewise linear ansatz (|1.7[) . u — ^2 rn k \x — x k \, 
is a weak solution to the linearly forced inviscid Burgers equation (|1.2[) if and 
only if 

n n 

Xk = y^ra* \xk - Xj\ , rhk = 2m fc J^m; sgn^ - x k ), (2.1) 

i=l i=l 

for k — 1, . . . , n. For this class of solutions, equation (|1.2|) takes the form 

u t + \{u 2 ) x = M 2 x - MM+, (2.2) 

where M = m k and M + — y^"_ 1 m k x k are constants of motion. 

Proof. This is a special case (all s k — 0) of Theorem [7] which is proved later. □ 

One can assume that all m k ^ 0, since it follows from (|2.ip that any vanishing 
mk remains identically zero. If we think of Xk and mk as positions and masses 
of particles on a line, then the total mass M and the center of mass M+ /M (if 
M 7^ 0) are conserved. Note that when M = we have the unforced Burgers 
equation. There are some additional constants of motion M2, . . . , M n that come 
together with Mi = M from the Lax pair presented in the next section, but we 
will not need them here [T]. 

The presence of absolute values and the sign function in (|2.1[) naturally 
divides the position space R ra into sectors. More precisely, to any permutation 
g = a\ era ... cr„ of the numbers {1,2, ... ,n} one can assign the sector X a — 
{(x%, X2, ■ ■ ■ , x n ) G R" I x ai < x a2 < ■ ■ ■ < x an }. We will concentrate on the 
sector X e corresponding to the identity permutation e =12. ..n, since there is 
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no loss of generality in assuming that the initial positions x k (0) are sorted in 
increasing order: 

X e = {{x 1 ,x 2 , ■ ■ ■ ,x n ) G R" \x l < x 2 < ■ ■ ■ < x n }. (2.3) 

For positions in X e the ODEs (|2.ip take the form 

n n 

ik = rnj(xk - Xj) sgn(fc — i), rhk = 2mfc rrij sga(i — k). (2.4) 
i=i i=i 

The following theorem solves this system completely. 

Theorem 2. Given any initial data {x k (0), rnk(0)}^ =1 (with the Xk{0) 's ordered 
or not), the solution of the ODEs (|2 .4[) is given by the formulas below, where 
M = ^2 m k and M+ — ^ m fe x fc as before, and where the empty sums Y]y and 
* n Fq an d Fn are t° be interpreted as zero (so that Fo(t) — e~ Mt and 
F n (t) = e Mt ). 



When M ^ the solution of (JUIJ) is 



\3<k 



m k (t) 



M 

TOfc(O) 



^(x fe (0)-^(0))m J (0) 



(2.5) 



F k ^{t)F k {ty 

for k — 1, . . . , n 7 where 



Mt I k \ — Mt I n 

• When M = i/ie solution of (12. 4| is 



= ^ (0) + * [ E M°) - ^ (°)) m i (°) - E (°) - x i (°)) m ^ (°) I > 



/.\ m fe(0) 



Fk-i{t)F k {t) ' 

(2.7) 

/or fc = 1, . . . , n, where 

(k n \ 

J=l j=fc+l / 
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• Letting Ik — Xk+\ — Xk for k = 1, . . . , n — 1, we have in both cases 



h{t) = h(0)F k (t). 



(2.9) 



The proof is presented at the end of this section. As an immediate corollary 
we obtain information about the original ODEs (|2.ip . 

Theorem 3. Given initial data {a^fe(O), mfc(0)}^ =1 to the ODEs (12. lj) such that 
Xi(0) < £2(0) < • • • < x n (0) (that is, with the positions in the sector X e ofH n ), 
the solution is given locally ( around t — 0) by the formulas of Theorem [Jl and 
this solution is valid as long as the positions Xk{t) remain in X e . 

A local solution that starts in X e hits the boundary of X e whenever Xk = 
Xk+i for at least one fc, an event which we refer to as a collision. It is clear from 
(|2.9[) that a collision occurs when some Fk becomes zero, at which time mk and 
rrik+i blow up. The local solution is valid up until the time of the first collision. 
In general a shock will then form, and the continuous ansatz (|1.7|) will not be 
able to describe the solution beyond the point of collision. We will return to 
this in the section about discontinuous solutions. 

If all mfe(0)'s have the same sign, then (|2.6p shows that there are no collisions, 
so the solution is global. In the case when all are positive, the asymptotic 
behaviour of this global solution as t — > +00 is that xi — > M + /M and mi — ► M, 
while Xk — ► +00 and mk for all k > 1. When the mfc(0)'s have mixed 
signs, collisions may or may not occur for t > 0. For example, in the case 
n = 2 a collision takes place when F\ (t) = (mi(0)e M * + m2(0)e~ M *)/M becomes 
zero, which happens when m2(0)/mi(0) < and t — (2M) _1 In m2(0)/mi(0)|. 
Consideration of cases shows that this value of t is positive iff mi(0) < < 
m 2 (0). 

The event when lk-i = h = is called a triple collision, since three particles 
come together at one point. The absence of triple collisions in the CH equation is 
a nontrivial result 8, 13J, but for the linearly forced Burgers equation it is much 
simpler. (Note, however, that triple collisions are possible for discontinuous 
piecewise linear solution; see the examples at the end of the paper.) 

Theorem 4. Collisions occuring in continuous piecewise linear solutions of the 
linearly forced Burgers equation (|1.2p cannot be triple collisions. 

Proof. A triple collision would occur if lk—i(to) = = £fc(to) for some to 7 which 
amounts to Fk-i(to) = = Fk(t ) by (|2.9[) . From the definition of Fk, it is 
obvious that this is impossible in the case M = 0, since we are assuming mk 7^ 0. 
In the case M ^ 0, it is also impossible, although less obvious; Fk(t ) = iff 

to = oif l°g — fnfif* " b ^°l nl and the quotient inside the logarithm is positive. 



So F k = F fe _! = iff = < 0, where A = E^ro^O) and 



rnk{A + mk + B) = 0. But this is ruled out by mk and A + mk + B = M both 



>j>k 3 




being nonzero. 



□ 
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We finish this section with the postponed proof of the main theorem. 



Proof of Theorem^ Assume to begin with that xi(0) < ••• < x n (0). Then 
(|2.4p is equivalent to (|2.1[) . and we can attack the problem by trying to find 
Xk(t) and TOfc(i) such that the corresponding piecewise linear u{x,t) given by 
(jl.7p satisfies the PDE (|2.2p . The x^'s divide the real line into n + 1 intervals 
which we number by k = 0, . . . , n. In each such interval u takes the form 
u(x,t) = a k (t)x + bk(t). Inserting this into (|2.2p yields a k + a\ = M 2 and 
bk + &feafe = —MM+, from which (in the case M^O) 

(t] = M M0) cosh(Mt) + Msinh(Mt) 
fcU a fe (0)sinh(Mt) + M cosh(Mt) { ' 1 

is found immediately, and by making an ansatz for b k with the same denominator 
as a k one also obtains 

b = a fc (0)M + (l - cosh(M^)) + M(b k (0) - M+ sinh(Mt)) 
fetl ^ a fc (0)sinh(Mt) + Mcosh(Mi) ' ^' ' 

Now Xfe(t) and mfe(t) are recovered from the relations mfc = ^(a^ — ak-i) and 
Xfc = — [bk — &jt-i)/(afe — ak-i). Because of the algebraic nature of the formulas 
thus obtained, they satisfy the ODEs (|2.4p identically, which shows that the 
assumption x\ < ... < x n is immaterial and can be removed. (This will be 
important later; see the comments after Theorem [5J) The simpler case M = 
(unforced Burgers) is entirely similar, except that 

afcW = M0)TT' 5fcW = M0)TT- (2 - 12) 

(The solution for M = can also be obtained by expanding e ±Mt = 1 ± Mt + 
0(M 2 ) in the solution for M ^ and letting M 0.) □ 



3 Inverse spectral construction of solutions 

The Lax pair 

-dl(t> = zm<f>, (3.1) 
(l>t=[z- 1 d* + c + u x -ud x ]<f> ) (3.2) 

with c an arbitrary constant, is compatible iff nit + m x u + 3mu x = and 
m x — u xxxi under the assumption of sufficient smoothness needed to justify the 
cross-differentiation. In particular, it is compatible if u evolves according to the 
derivative Burgers equation (jl.ip . which can be written as m t + m x u + imu x — 
with To obtain the linearly forced Burgers equation (|1.2p from 

equation (jl.ip the rule (u 2 ) x = 2uu x is used. It is not obvious if all these 
formal calculations have any relevance to weak solutions, where the smoothness 
assumptions may be violated. To investigate this, let us say that (|3.ip and (|3.2p 
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constitute a weak Lax pair if they are satisfied in the weak sense discussed in the 
introduction (thus (j), like u, is a 2?'(R)-valued function of t, and the equations 
hold in the space of distributions X>'(R)). Solutions u of the form (|1.7[) . u = 
^2 rrik \x — xu\, do admit a weak Lax pair with m = u xx — 2^?, mfc&cj., 
and <j> is in this case a continuous function (in fact, it is piecewise a quadratic 
polynomial in x with i-dependent coefficients) . The product m<f> in (|3.1|) is well- 
defined since the distribution m can be multiplied by the continuous function <f>. 
We hope to treat weak Lax pairs in more depth in future papers. Here we 
just state a theorem which can be verified by careful use of the calculus of 
distributions. 

Theorem 5. The following are equivalent conditions on a junction u of the 
form ljl.7j) . u — Y, m k \x - Xk\: 

1. u is a weak solution to the linearly forced Burgers equation (|1.2p . and 
{xk , rrik } satisfy equations (|2.1[) . 

2. u has a weak Lax pair (|3.1|) , (|3.2|) . 

When u = ^2 m k \ x ~ x k\i a solution to equation (|3.1[) with the asymptotic 
condition (f>(x,t;z) — 1 for x < X\(t) will be consistent with the time evolu- 
tion given by (|3.2p provided that we choose the constant c = — M . Such a 
solution evaluated at x > x n (t) will take the form (f>(x,t;z) = A{t;z)\(x — 
x n ) 2 + B(t; z)(x — x n ) + C(t; z), where all three coefficients are polynomials in 
z, which, by equation ([3T2]) . satisfy A = 0, B = MB, and C = 4 + 2M C (see 
[lj. Thus it is consistent with equations (|3.1[) and (13. 2p to impose the condition 
A(t] z) — 0, which together with <f> = 1 for a; < x\ amounts to the boundary 
conditions 4> x {~ oo) = <fi xx (—oo) = <f) xx {oo) = 0. With these boundary condi- 
tions in place, the problem of solving the ODEs (|2.1[) becomes an isospectral 
deformation problem which can be solved if one knows how to solve the inverse 
problem for equation ()3.1|) . This is exactly the inverse problem that was studied 
in [1] under the additional assumption that all mfc(0) > 0. We now give a brief 
summary of results from that paper. 

Theorem 6. The "Neumann-like discrete cubic string" boundary value problem 

-d x 4> = zrrxj), <f> x (-oo) = cj> xx (-oo) = c/) xx (oo) = 0, 

where m — 2J^ =1 mkS Xk with all > 0, has a spectrum of the form {0 = 
zq < z\ < Z2 < ■ ■ ■ < z n -i}. There is a one-to-one (up to translations of m 
along the x axis) and onto spectral map m i— > {M, ji}, where M — > 

and [i is a measure of the form /i = X)j=i bj8 Zj , with bj > for j = 1, . . . , n — 1 
(see details in J^j). The inverse problem of recovering the discrete measure m 
from {M, ji\ has the explicit solution 

_ CkVk _ _ , _ 2Ak „s 

^n — k — . A 3 %n — k+1 %n — k — <"n — k — • v^-^J 

ZAk+lAk V k 

in terms of determinants of bimoment matrices constructed out of the measure 
[i and the constant M (see below). 



8 



We recall the following definitions from pQ. Given a measure /i, let 



pj = / z 3 dfi(z), 



z + to 



dfj,(z)dfj,(w). 



(3.4) 



Let .Ao = Bo = Co = V a = 1, Ai = ho + 517 ■ D'i = Po, an d for other values of 
k let 



B k 



I()Q 

ho 



A h 



hi 
hi 



ho + 2M 

ho 

1 20 



hi 
hi 
hi 



h-1,0 h-1,1 

ho hi 

ho hi 

ho hi 



Ik-1,0 Ik-1,1 

h,k-i 
h,k-i 

Ik-l,k-l 

h, k -i 

h,k-l 



h.k-1 



h,k-l 
h,k-l 
h,k-l 

Ik-X,k-l 

hi I12 
hi I22 

hi Ik2 

Po 1 10 
Pi ho 

Pk-1 ho 



I Ik 
I 2k 

Ikk 

h,k~2 
h,k-2 



'fe,fe-2 



(3.5) 



In all these cases, the index k agrees with the size k x k of the determinant. 
Note that A k = B k + sifCfc-i for fc > 1. 

Let us analyze the formula (|3.3p for h in order to compare it with (|2.9[) 
obtained earlier. First, (|3.2p implies that the linearly forced Burgers equation 
induces a very simple evolution of the measure /i, namely /i(z; t) = e Mt [i{z\ 0). 
Because of this it is easy to factor out the time dependence from all the deter- 
minants involved in (|3.3[) . This elementary exercise leads to h(t) — lk{0)h\{t), 
where 

_ B ,- 1 .(o)^ + Aft.- t -(Q) : -'". M 



Bn-k(0) + 3^C n _ fc _i(Q) 



This is in full agreement with 
checked in a similar way. 



and 



The formula for m k can be 



4 Discontinuous piecewise linear solutions 

Theorem 7. The discontinuous piecewise linear ansatz (|1.9p . u — XX m & \ x ~ x k\~ 
Sk sgn(a: — Xk)), is a weak solution of the linearly forced inviscid Burgers equation 
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(11. 2[) zf and onfr/ if 



n 

±k = }Xmi \xk -Xi\ + s t sgn(xi - x k )), 

1=1 (4-1) 
m k = 2m k ^ m t sgn(xi - x k ), s k = s k ^ m i sg n ( x i ~x k ), 

i=l i=l 



for k — 1, . . . , n. For this class of solutions, equation (jl.2|) takes the form 

u t + \(u 2 ) x = M 2 x - M(M+ + S), (4.2) 

with M = ^2 m k and M+ = ^2m k x k as before, and with S = ^2s k . The 
quantities M and M + + S are constants of motion, and so is s\jm k for k = 
1, .... n (provided that m k ^ 0). 

Proof. We will repeatedly use the following distributional formula valid for 
an arbitrary piecewise differentiable function / with points of discontinuity at 
X\ , Xi , . . . , x n 

'■ fx = {fx} + J2 k =i[f]k5 Xk , where {f x } means the ordinary deriva- 
tive taken away from discontinuities and [f] k — f(x^) — f{xZ) denotes the jump 
at x k . Moreover, u t = J2 k (m k \x - x k \ - (m k x k + s' k ) sgn(x - x k ) + 2s k x k 5 Xk ). 
Now the left-hand side of (|1.2|l . Ut + ^{u 2 ) x , must be a function since the right- 
hand side is a function; hence all Dirac deltas must cancel out. Similarly, there 
must be no Dirac deltas in the first or second x derivatives of Ut + ^(u 2 ) x , These 
conditions give, in turn, 



= 2s k x k + \[u 2 ] k , = -2(m k x k + s k ) + \[{{u 2 ) x }] k , 
Q = 2m k + \[{{u 2 ) xx }] k . 



(4.3) 



An elementary computation of jumps for the case of piecewise continuous func- 
tions now produces (|4.ip . The coefficients of the forcing term in the PDE are 
identified from the smooth part of the term ^(u 2 ) x , while the constants of mo- 
tion follow from (jlTTj) . □ 

Weak solutions to an initial value problem are usually not unique unless 
the PDE is supplemented with a so-called entropy condition that picks out the 
"physical" solution. In the case of the Burgers equation this condition requires 
u to jump down, not up, at discontinuities. This is satisfied by the ansatz (|1.9| 
if all shock strenght s k are nonnegative, so we will assume s k > from now on. 

When considering the initial value problem for the ODEs (14.11) we can as- 
sume without loss of generality that cci(O) < X2{0) < ... < x„(0). Thus on a 
sufficiently small time interval we will still have xi(t) < x 2 {t) < . . . < x n (t); in 
other words, the positions stay in the sector X e (see (|2.3p ). In X e the equations 
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(|4.ip can be written as 



n 



x k = }Xmj sgn(fc - i)(xk - Xi) - s, sgn(fc - i)), 

1=1 (4.4) 

n n V / 

rfife = 2m fe m, sgn(i - /c), s^st^m, sgn(i - k), 

i=l i=l 

for k = l,...,n. These equations can be solved explicitly, since the simple 
change of variables in the following theorem reduces them to the ODEs already 
solved in Theorem [5] 

Theorem 8. // {xk, rrik, satisfy (|4.4p . if all m,k(0) ^0, and if 

Vh = %k + Sfc/mfc, (4.5) 
then {j/fe, rrik}k =1 satisfy (|2.4[) (with yk taking the place of Xk everywhere). 
Proof. Straightforward calculation. □ 

Note that the initital values ?/fc(0) will not necessarily be distinct or sorted in 
increasing order even though the cc/c(0)'s are, but this does not matter since the 
solution formulas of Theorem[2]are valid for any initial conditions. So Theorem^ 
gives us yk{t) and TOfc(i) (note that ^ mkyk = XX TO fc x fc + s fc) = M+ + S replaces 
M + in the solution formula (|2.5p ), and we can then recover Sk(t) from the fact 
that s\jm,k is constant for each k\ this gives Sk(t) = Sk(0)/ \J Fk-i(t)Fk(t), and 
allows us to also recover Xk(t) — yk{t) — Sk(t) / rrik(t) . This solution {xk, mk, Sk} 
to (14. 4p is also the solution to (|4.1|) , at least locally in some time interval around 
t = (so that the x/c's remain in the sector X e ). 

For illustration, here is the general solution with shocks in the case n = 2, 
when M j£ 0, mi(0) ^ 0, m 2 (0) ^ 0: 

, . "ii(0) . . si(0) 

mi t = — —7— — ——— , s x it) = , - 

m w 2(0) m * 2 (0) 

m 2 (t) = 777-777777, s 2 {t) - 



Fi(t)F 3 (t)' y/F^F^ty 
M+ + S - Km 2 (0)e~ Mt si(Q) 



M mi(0) ^ 46 ) 

M + + 5 + ifmi(0)e Mt s 2 (Q) 



* i(i > = - M 

*-„<«)-„<<»+ * m s,(0) 



A// 



m 2 (0) mi(0) 



In the continuous case (|2.1[) we assumed all mk 7^ 0, but for (|4. 1[) it does 
make sense to have rrifc = provided that the corresponding Sk is nonzero. If 



11 



-«(*) o 



?(*) 



X 



X 








Figure 1: Left /middle: Wave profile u(x,t) as given by (|4.7p at two different 
times i < t co u, with decreasing towards zero at a constant rate. Right: 
Stationary profile after collision (i > £ co n). 

mfe(O) = 0, then clearly mu{t) — for all t, and the above solution procedure 
does not work. But this is easily fixed: just write down the general solution 
obtained for a nonzero initial value mfe(0) = a, and let a — > there. 

We will finish with a few examples that show how to deal with the solution 
when it hits the boundary of the sector X e . 

Example. A particular antisymmetric solution of (|4.1jl with n = 3 is given 
by —x\ = xs = £ > 0, X2 = 0, —mi = m 3 = fi > 0, m 2 = 0, si = s 3 = 0, 



s 2 = a > 0, where £(t) = ^(0)F(t) - cr(0)t, /i(t) = ^(Q)/F(t), a{t) = a(0)/F(t), 



with F(t) = 1 — 2/i(0)t. (These formulas are obtained either by reducing (|4.ip 
to ODEs for £, /x, cr and solving them directly; or by assuming TO2(0) = a / 0, 
changing variables to y\ = xx, j/2 = X2 + S2/rri2, 2/3 = £3, writing down the 
general solution using Theorems [8] and [21 and letting a — > 0; or simply by 
noting that M = so that we are dealing with the unforced Burgers equation 
whose solution can be found in the textbook way using characteristics.) Since 
M + + S = 2/i£ + a is constant in time, the wave profile (see Figure [T]) is 



If cr(0) = then this is a shockless solution (with n = 2 really, since there is 
neither mass nor shock at the site X2 = 0). It is defined until = £(0)F(t) 



while mi and 7713 blow up to —00 and +00, respectively. However, u remains 
bounded, and tends to a shock profile: u(x, t) — > —2^(0)^(0) sgn(x) as t / t co \\. 
This illustrates that shocks can form naturally even if they are not present 
in the initial wave profile. The profile will be stationary after the collision, 
because its continued evolution is given by the n — 1 case of (|4. 1[) (x\ = mi, 
rhx = si = 0) with x\ = 0, mi = 0, Si = 2/j(0)£(0). Consequently, u(x,t) = 
-2/x(0)£(0) sgn(x) for all t > t coi i. 



u(x, t) 



-flit) \x + I + n(t) \x - £(t)\ ~ <?(t) sgn(x) 

'2/i(0)€(0)+(r(0), x<-{(t), 

-2fi(t)x + a(t), -£(*)< a; <0, 

< 0, a; = 0, 

-2fi(t)x - a(t), 0<x<£(t), 

k -(2/i(Q)^(0)+tr(0)), f(t)<a;. 



(4.7) 



becomes 




Then a;i and X3 collide at a; = 
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xi(0) = -2 a; 2 (0) = x 3 (0) = 1 x 

Figure 2: Solid: Continuous initial wave profile u(x, 0). Dashed/dotted: u(x,t) 
at times t = In | and i = In | , respectively. 

If <t(0) > there is a shock waiting at the origin between the two approaching 
particles (as in Figure [TJ. The solution hits the boundary of the sector X e when 

£(t) becomes zero at time t co ii = (2/i(0)+cr(0)/£(0)) . Then xi = x 2 = 2:3 = 0, 
which illustrates that triple collisions may occur when shocks are present. Since 
the collision occurs earlier than in the shockless case, F(t) has not yet reached 
zero at the time of collision; hence mi and 7713 do not blow up in this case. Again, 
u tends to a stationary shock profile: u(x,t) = — (2/x(0)£(0) + <r(0)) sgn(x) for 
all t > i C oii- 

Example. Consider now the shockless ODEs (|2.1I) with n — 3 and initial data 
mi(0) = |, 7712(0) = —1 and m 3 (0) = |, so that M = 1. We assume xi(0) < 
^2(0) < £3(0) but leave them otherwise unspecified. Since u — ±(Mx — M + ) 
as x — > ±00, and since the slope u x jumps by 2rrik at each Xk, the initial 
profile u{x, 0) consists of line segments with slope —1, -|, — | and 1, joined 
at the points [x^, u(xk, 0)). Figure [2] illustrates this for the particular values 
xi(0) = -2, x 2 (0) = 0, x 3 (0) = 1. Note that if the lines u = ±(Mx - M+) to 
the left and to the right are continued, they intersect on the x axis at the center 
of mass x — M + /M (= in the figure), which is a constant of motion. 

Recall that h = X2 ~ x\ and I2 = £3 — x 2 . From (|2.6|) we obtain Fo(t) = e - *, 
Fi(t) = fe* + \e-\ F 2 {t) = -|e* + |e -t , and F 3 (t) = e*. Equations dlHJ and 
(EU) givex 1 (t)=a; 1 (0) + (l-e- t )(iii(0) + f/2(0)), i 2 (t) = n (t) + Zi(0)Fi (t), 
and xs(t) = x 2 (t)+/ 2 (0)F2(7j). There is a collision between £2 and X3 when i*2(f) 
becomes zero, which happens at time t = t co n = In 2 when e = 2. At that time 
we have F = |, F x = |, F 2 = 0, F 3 = 2, hence by (EJ3 mi = m 1 {0)/F o F 1 = §, 
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zi(toii) = -1 x 2 (t co n) = x 3 (t co n) = 1 x 

Figure 3: Solid: Discontinuous wave profile u(x,t) formed at the instant of 
collision t = t C oii = In 2. Dashed: u(x,t) at time t = t co n + \ . Dotted: No more 
collisions occur, and u(x, t) — > |x| as t — > +oo. 

77i2 = — oo, 7713 = +00. As for the wave profile u, we have 

u(xi(t),t) — mil\ + 1713(1% + h) = (M - mi) h + TI3Z2 

= (Fi - mi(0)/F ) h(0) + m 3 (0)Z 2 (0)/F 3 (4.8) 
-» I ii(0) + f i 2 (0), as i co ii, 

and 

It(x2(t), t) - u(x 3 (i), t) = (mill + m 3 Z 2 ) — (mi(!i + Z 2 ) + TO2Z2) 
= (m 3 — mi — m 2 ) Z2 

= (m 3 (0)/F 3 - m 1 (0)F 2 /F F 1 - m 2 (0)/F 1 ) l 2 (0) 
->§Z 2 (0), ast/'t co ii. 

(4.9) 

Thus the limiting wave profile at i = t co n consists of a line segment with slope 
— 1, joined to a line segment with slope —1 + 2- | = |ata: = a;i(t C oll) = 
»i(0) + |(^i(0) + 3/2(0)) and height u= \ h(0) + § Z 2 (0); the profile jumps 
down by |Z 2 (0) at x = x 2 (t C oii) = X3(t ou) = £i(*coii) + f ^i(O), and continues 
from there with slope 1. See Figure [3] 

The continued evolution of the profile for t > t co \\ is illustrated in Figure[3J it 
is given by the shock ODEs (|4. 1[) with n = 2, using a new set of variables whose 
initial values at t = t co \\ are x\ — xi(t co \\), x 2 — £ 2 (i C oii), Thi — §, TO 2 = §) 
Si = 0, and s 2 = §^2(0). In terms of the new time variable r — t — t co u > one 
finds from the general solution (|4. 6[) that, for example, 

x a (r)-£i(r) - f 5 2 (Q) - £ x (0) + 4^) Fi(r) -4^^( T ^W. ( 4 - 10 ) 
V m 2 (0)/ m 2 (0) v 
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where Fi(t) = |e T + T and ^(t) = e T . Writing this expression as x 2 — x,\ = 

(A+B)F 1 -By / F 1 F 2 , we see that it is zero if F± (r) = 0, which can never happen, 
or if (A + B) 2 Fi = B 2 F 2 , which is the same as e~ 2T = 9((A + B)/B) 2 - 8 that 
can't happen either since the right-hand side is > 1 and the left-hand side is 
< 1 for t > 0. The conclusion is that, in this example, x 2 (t) — x\(t) remains 
positive for all r > 0, so there are no more collisions. Instead, as r (or t) — > +oo, 
we have x\ — > 0, x 2 — > +oo, mi — ► M, rh 2 — ► 0, and s 2 — * 0. Thus, u(x, t) 
approaches the limiting wave profile u(x, +oo) = \x\. 
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